Magneto-orbital helices and multiferroicity 
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Orbitally ordered states evolving coincidentally with structural distortions and 
magnetic ordering provide a novel route to controlling electronic and magnetic prop- 
erties of materials through external fields. We report an unprecedented magneto- 
orbital helix in CaMn70i2, found to give rise to the largest magnetically induced 
ferroelectric polarisation measured to date. Characterisation of the structural mod- 
ulation using x-ray diffraction, and analysis of magnetic exchange shows that orbital 
order is crucial in stabilising a chiral magnetic structure, thus allowing for electric 
polarisation. Additionally, the presence of a global structural rotation enables the 
coupling between this polarisation and magnetic helicity required for multiferrocitiy. 
These novel principles open up the possibility of discovering new high-temperature 
multiferroics. 

Orbital ordering phenomena provide one of the most striking manifestations of electronic corre- 
lations in a variety of materials, most notably in 3d transition metal oxides pp. In these systems, 
below a certain temperature Too, a distinct anisotropic pattern of orbital occupation emerges 
from a nearly isotropic high-temperature state. Since this state is strongly connected with both 
structural distortions (the cooperative Jahn Teller effect) and with magnetism, orbital ordering 
has been regarded for many years as a promising ingredient to allow the cross-coupling of electric 
properties. For example, colossal magnetoresistance takes place because an orbitally ordered, in- 
sulating and antiferromagnetic state competes with a ferromagnetic isotropic metallic state, and 
its onset can be controlled by the application of an external magnetic field [2]. Recently, it has 
been proposed that some charge and orbitally ordered states could be electrically polar, so that 
orbital ordering could promote high-temperature, magnetically controllable ferroelectricity [3]. 
Here, we show that a unique form of incommensurate orbital ordering in the complex manganese 
oxide CaMn70i2, appearing below 250 K, is responsible for the largest observed magnetically in- 
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duced ferroelectric polarisation. Unlike all other known systems, orbital occupation in CaMn70i2 
does not select fixed preferential directions, but rotates around the 3-fold crystallographic axis, 
preserving the overall point-group symmetry. The resulting pattern of magnetic exchange inter- 
actions produces a magnetic helix having twice the period of the orbital modulation below 90 
K. This magnetic texture breaks inversion symmetry and couples to the crystal structure via the 
relativistic antisymmetric exchange to yield a large electrical polarisation. 

Figure la shows the room temperature crystal structure of CaMnyO^. It can be described 
as a perovskite superstructure, with 3/4 of the A sites occupied by Mn 3+ (hereafter labelled 
Mnl), and the remaining occupied by Ca. A trigonal structure (space group R3) develops on 
cooling below ~ 400 K, and is best described by considering the pseudo-cubic B sites, which form 
three inequivalent triangular layers stacked along the hexagonal c axis (shown in Figure la, and 
coloured green, blue and red throughout this report). Within each layer, 1/4 of the B sites retain 
3 symmetry (Mn3, black-outlined circles in Figure la), whereas the other 3/4 have 1 symmetry 
(Mn2). Charge ordering develops gradually on further cooling, with Mn2 and Mn3 tending towards 
Mn 3+ and Mn 4+ , respectively. An incommensurate structural modulation along the hexagonal c 
axis appears below T oo =250K, with a propagation vector q c = (0, 0, 2.077±0.001) at 150 K, 
producing a continuous variation in Mn-0 bond lengths [7J. 

The magnetic structure below Tn = 90 K was determined using neutron powder diffraction [I] , 
and is shown in Figure lb. Magnetic anisotropy results in the spins lying within the hexagonal ab 
plane. Within each of the triangular layers, Mn2 spins are parallel to each other, whilst Mn3 spins 
are approximately antiparallel. Adjacent layers are rotated by 124 ° in the ab plane with respect 
to each other, leading to an incommensurate magnetic modulation with propagation vector q m 
= (0, 0, 1.037) (Supplementary Section SI). As observed in many magnetically ordered crystals 
(see, for example [3]), below T N , q c = 2q m . In the present case, the structural modulation is also 
present in the paramagnetic phase, and the magnetic modulation "locks in" to it as it develops at 
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We now focus on the incommensurate structural modulation appearing below Too, and demon- 
strate that this corresponds to the onset of a novel form of orbital ordering. We performed Rietveld 
refinement using JANA2006 [6J, of single crystal diffraction data measured at 150 K (Supplemen- 
tary Section SII). Consistent with previous reports [7], we found that the modulated crystal 
structure is described by the centrosymmetric four-dimensional space group R3 (007)0 with prop- 
agation vector q c = (0, 0, 2.077). In Fig 2a, we plot the Mn-0 bond lengths for adjacent unit 
cells along the c axis around Mn2 (the variation in the Mnl-0 and Mn3-0 bond lengths are small 
and dominated by movement of the oxygen anions (Supplementary Section SIII)). This structural 
modulation has strong implications for the nature of orbital occupation as follows. We define 
orthogonal x, y and z bond directions, schematically shown in Figure lc. In the unmodulated, 
average structure, the Mn 3+ 06 octahedra around Mn2 are compressed, giving two short and four 
long bonds, implying an unusual occupation of the x 2 — y 2 orbital. In reality, the structural 
modulation leads to an incommensurate rotation of a single pair of long bonds in the xy plane, 
consistent with a periodic occupation of the 3x 2 — r 2 and 3y 2 — r 2 orbitals, shown at the side of 
Figure 2. We can further characterise this orbital rotation by employing a formalism [S] [DJ, first 
proposed by Goodenough, which quantifies orbital occupation according to a mixing angle, tan# 
= y/3(x — z)/(2y — x — z). Following the conventions in [9], we need only to consider the angular 
sector 27r/3 < 9 < An/3 (inset to Figure 2). As shown in Fig 2b, the orbital occupation rotates 
periodically between the 3x 2 — r 2 and 3y 2 — r 2 orbitals, passing through the x 2 — y 2 only at the 
nodal points. Figure 2c shows the Mn2 charge modulation accompanying the orbital modulation, 
calculated using the Bond Valence Sum method [ID] . Valence fluctuations are surprisingly small 
( 1 % — Supplementary Section SIV), indicating that the modulation affects predominantly the 
orbital occupation. 

Although orbital ordering is a widespread phenomenon in manganese oxides jTl] [12], an in- 
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commensurate modulation where the orbital occupation rotates perpendicular to the propagation 
direction is unprecedented. The resulting orbital arrangement has a profound effect on magnetic 
ordering below TV = 90K. This can be seen by applying the Goodenough-Kanamori-Anderson 
(GKA) rules [13J [H] [15] [16], which predict the nature of exchange interactions. In the absence of 
the structural modulation, the pseudo-perovskite structure is heavily distorted, and antiferromag- 
netic interactions are weak due to the large departure of the Mn-O-Mn bond angles from 180°. 
The alternation of filled and empty x 2 — y 2 orbitals on Mn2 and Mn3 promotes ferromagnetic 
exchange, and a net ferromagnetic structure is predicted. The situation changes drastically in 
the presence of the orbital modulation, as the magnetic interaction between adjacent A, B and 
C layers is predominantly mediated by the orbitals that are most affected by the modulation. 
This is shown schematically in Fig 3a, which depicts a region of the crystal structure where the 
modulation amplitude is largest. Mn3 in the B layer is connected to the Mn2 sites in the A layer 
by filled 3x 2 — r 2 orbitals, yielding a ferromagnetic interaction J 1; and to the Mn2 sites in the C 
layer by empty orbitals, yielding an antiferromagnetic interaction J 2 . This produces an effective 
ferromagnetic interaction between spins in each of the A, B and C layers, as all Ji(J2)'s are equiv- 
alent by symmetry (shown in Figure 3b). Across layers there is an effective next-nearest neighbour 
antiferromagnetic interaction J12, which is strongest at the "crest" of the modulation, but remains 
antiferromagnetic everywhere except possibly at the nodal point. The overall nearest-neighbour 
interaction between layers, J3, results from the competition between full-empty (x-z) orbitals 
and empty-empty (y-z) orbitals, which are ferro- and antiferromagnetic, respectively. Although 
the sign of J3 cannot be predicted a priori, it is experimentally found to be positive (antiferro- 
magnetic). The magnetic helix therefore results from the competition between antiferromagnetic 
nearest-neighbour and next-nearest neighbour interactions. We can test this hypothesis with a 
very simple model, which captures the effect of the average interactions Ji, J2 and J3, assuming a 
magnetic propagation vector (0,0, q 2 ) and identical spins S on Mn2 and Mn3 sites. The magnetic 
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Heisenberg exchange energy per Mn3 site can be written as: 

E = S 2 [3JiCos0 + 3J 2 cos(>-2^) + 6J 3 cos^] (1) 

where tp = 2/37rg z and are the angles between Mn2-Mn2 and Mn3-Mn2 spins in adjacent 
layers, respectively. Minimising the energy in eq. [T] yields the phase diagram shown in Figure 
3c. The wide incommensurate phase field is bound on either side by two distinct commensurate 
collinear phases, both with q z = 1.5. The white line indicates the range of parameters for which 
the most stable propagation vector is that determined experimentally. The red dot corresponds 
precisely to the experimental magnetic structure (0 = 31°), and is found for J2/Ji= -0.85 and 
J3/J1—O.63, in good agreement with our qualitative considerations based on the GKA rules. The 
coupling between Mnl and adjacent Mn2 sites is expected (and observed) to be ferromagnetic 
(Figure 3a) since the sites are linked via oxygen anions with highly-distorted bonds at less than 
the GKA critical angle. The most important ingredient not included in our simple model is the q c 
= 2q m lock-in term, which intertwines the magnetic and orbital modulation to form a magneto- 
orbital helix. Phenomenologically, different forms of the lock-in term are allowed by symmetry, with 
the most likely candidate involving a third-harmonic " squaring up" of the magnetic modulation. 

Below T^, CaMn 7 0i2 acquires the largest measured magnetically-induced electric polarisation 
along the c axis [I], which can be explained by the existence of the magneto-orbital helix. As 
CaMn 7 0i2 has ferroaxial point symmetry 3, the electrical polarisation can arise through coupling 
between the magneto-orbital helix and a global rotation between different elements of the crystal 
structure, as shown in Figure 4a: the 6 oxygen atoms surrounding an individual Mn3 are all rotated 
in the same clockwise direction. This is described by the phenomeno logical invariant ao~A • P, 
where a is a coupling constant, a is the magnetic helicity, P is the electrical polarisation and A 
is an axial vector representing the structural rotation. The magnetic helicity is not invariant by 
a global rotation of the spins, so the coupling can only be of relativistic nature. The relevant 
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mechanism is the antisymmetric Dzyaloshinskii-Moriya (DM) exchange, where energy can be 
gained by distorting the crystal structure and/or the electronic density in the presence of non- 
collinear spin configurations. Although both structural and electronic distortions are allowed 
in the present case, we will employ the former to illustrate how this mechanism can produce 
electrical polarisation in CaMnyO^. When two metal sites carrying non-collinear spins are joined 
by common ligand atoms (in this case Mn-O-Mn) , energy can be gained by displacing the ligand 
through a vector u so that AE = Au ■ [ri 2 X (Si x S2)] = AD • (Si x S2), where Si and S2 
are the spins on the two sites, r 12 is the position vector connecting them and A is a coupling 
constant. D = u x ri2 is the familiar DM vector |TT] [18] - Energy is therefore minimised by a 
pattern of local ligand (oxygen) displacements u associated with pairs of spins and either parallel 
or antiparallel to the vectors ri 2 x (Si x S2), depending on the sign of A. Fig 4b shows the pattern 
of ri2 x (Si x S2) vectors (black arrows) for clusters of Mn2 around a single Mn3. Assuming 
A < 0, the DM interaction will favour oxygen displacements parallel to these arrows, making the 
Mn 3+ -0-Mn 4+ bonds between the red and the blue layer flatter (i.e., the bond angles closer to 
180°) and Mn 3+ -0-Mn 4+ bonds between the red and the green layer more acute (i.e., the bond 
angles farther from 180°). It is easy to see that this pattern of displacements will result in the 
central Mn 4+ moving upwards (yellow arrow), creating a local polarisation along the positive c 
direction (identical considerations for A > lead to a polarisation along the negative c direction). 
Since the magnetic structure is globally chiral, these local polarisations add up to give a net 
electrical polarisation. It is important to stress that such a simple relation between local u and 
global P vectors is only possible in the presence of a global structural rotation, which is allowed 
in point group 3, as no net polarisation would arise for 180° Mn 3+ -0-Mn 4+ bonds. Ferroaxiality 
is thus an essential ingredient for the development of electrical polarisation in CaMnyO]^. The 
large value of P is readily explained by considering the angle of rotation <p between the spins 
Si and S2, since P should be roughly proportional to sin0. The average angle of rotation is 
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(0) = 124/2° = 62° whereas for the prototypical magnetic multiferroic TbMn0 3 , <p = 200°, 
so that sin(0) in CaMn70i2 is 2.5 times larger than in TbMnOs, in relatively good agreement 
with that found experimentally (P four times larger), considering the differences in geometry and 
coupling. The small deviation from commensurability is not an important factor in determining 
the electrical properties of CaMnyO^: P would be almost as large if the magnetic propagation 
vector was exactly (0,0,1), since what matters is the large rotation angle between spins in adjacent 
A, B and C layers. 

To conclude, we have shown that the magnetic and structural modulations in CaMn 7 0i2 are 
intertwined to form an incommensurate magneto- orbital helix . This unprecedented texture rep- 
resents a novel manifestation of the deep connection between orbital physics and magnetism in 
transition-metal oxides. The incommensurate structural modulation is crucial in stabilising a 
chiral magnetic structure, which breaks inversion symmetry, producing the largest magnetically- 
induced electrical polarisation ever observed in any system. We have also shown how the pres- 
ence of a global structural rotation (ferroaxiality) is essential in promoting the coupling between 
magneto-orbital helicity and ferroelectricity. These observations open up the possibility of finding 
a new class of multiferroic materials in which giant ferroelectric polarisations may arise, under- 
pinning the development of multiferroic technology. 
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Figure 1 The crystal structure, magnetic structure and associated orbital orientations 
of CaMn 7 Oi2. (a) (left panel) The room temperature crystal structure in the hexagonal setting, 
highlighting the three triangular layers (A=blue, B=red and C=green) stacked along the hexag- 
onal c axis (cubic [111] direction). The structure can also be described within the pseudocubic 
perovskite setting (right panel), with Mn3 sites (black-outlines spheres, with nominal valence 4+) 
occupying the apex sites along the [111] diagonal, (b) The magnetic structure at 90 K. The tri- 
angular layers are again shown as blue, red and green, with a 124° spin rotation between adjacent 
layers, (c) A schematic showing the B sites in the pseudocubic setting, defining the x, y and z 
bond directions which are used to describe orbital orientation. 
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Figure 2 The incommensurate structural modulation coupled to orbital rotation and 
valence fluctuations (a) Variation in Mn2-0 bond lengths along x, y and z (coloured according 
to the schematic in Figure lc) as a function of R z , plotted across 19 unit cells along the c axis, 
(b) Orbital mixing angle 9 for the Mn2 octahedra as a function of R z . 6 correlates directly with 
orbital occupation, as shown by the circular inset. The orbital rotation is shown, along with the 
coupled helical spin rotation, (c) The Mn2 charge modulation, which accompanies the structural 
modulation along the c axis. The incommensurate orbital modulation for half a period of the 
structural modulation is depicted at the side of the figure, with positions 1,2 and 3 corresponding 
to the labelled points in (a) . 
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Figure 3 The effect of orbital ordering on magnetic exchange interactions, promot- 
ing a chiral magnetic structure (a) Three pseudocubic unit cells illustrating the coupling 
between the incommensurate orbital modulation and the magnetic structure. The occupied Mn2 
(Mn 3+ ) orbitals are shown in grey, giving rise to ferromagnetic exchange interactions (Ji), shown 
by a dashed line, and antiferromagnetic exchange interactions (J2), shown by a dotted line. The 
nearest-neighbour and effective next-nearest neighbour antiferromagnetic exchanges between lay- 
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ers are labelled J 3 and J 12 , respectively. The ferromagnetic exchange along the Mnl — Mn2 chains 
is shown by a solid black line, (b) Figure 3a viewed down the c axis, showing the symmetry 
equivalent exchange interactions, (c) Phase diagram of the CaMn 7 0i 2 magnetic structures, ob- 
tained by minimising the classical Heisenberg exchange energy in eq. 1 as a function of the angles 
ip and for different exchange constants Ji, J 2 and J3. Spins on Mn2 and Mn3 sites were taken 
to be equal (different spins would be equivalent to renormalising the exchange constants). CM1 
and CM2 are two commensurate collinear phases, both with propagation vector (0, 0, 1.5), in 
which the Mn2 and Mn3 in each layer are antiparallel (for CM1) and parallel (for CM2). ICM is 
an incommensurate phase with rotating spins (magnetic helix). The shading indicates the value 
of the ^-component of the propagation vector. The white line indicates the locus of points with 
the experimental propagation vector (0, 0, 1.037). The red dot corresponds to the experimental 
magnetic structure (0=31°). 



14 




Figure 4 The coupling between magnetic helicity and ferroelectricity, leading to a 
giant polarisation, (a) The global structural rotation described by the axial vector A along 
the c axis is evident in the pattern of oxygen atoms (red spheres), which couples to the clockwise 
(anticlockwise) spin rotation between Mn 4+ and green (blue) Mn 3+ layers, (b) A perspective 
view on the Mn 3+ — Mn 4+ — Mn 3+ configuration. The inverse Dzyaloshinskii-Moriya interaction 
favours a displacement of the oxygen anions in the direction of the black arrows. This results in 
an additional displacement of the Mn 4+ cation, indicated by the small yellow arrow, and thus a 
net polarisation along the positive c axis, shown by the large yellow arrow. 



15 



